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ON . We theoretically investigate resonant tunneling through a linear array of quantum dots with 

subsequent tunnel coupling. We consider two limiting cases: (i) strong Coulomb blockade, where 
only one extra electron can be present in the array (ii) limit of almost non-interacting electrons. We 
develop a density matrix description that incorporates the coupling of the dots to reservoirs. We 
analyze in detail the dependence of the stationary current on the electron energies, tunnel matrix 
elements and rates, and on the number of dots. We describe interaction and localization effects on 
the resonant current. We analyze the applicability of the approximation of independent conduction 
channels. We find that this approximation is not valid when at least one of the tunnel rates to the 
leads is comparable to the energy splitting of the states in the array. In this case the interference of 
' conduction processes through different channels suppresses the current. 
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<D I I. INTRODUCTION 

In recent years arrays of quantum dots have received an increasing amount of interest. With|-,the progress of 
fabrication techniques quantum dot arrays are coming within the reach of experimental investigations If the electron 
levels in the individual dots are aligned we encounter here a situation of resonant tunneling. In this regime, the 
1 , ■ transport in the array becomes sensitive to precise matching of the electron levels in the dots that can be controlled 
by external gates. This opens u» new possibilities to control the transport and perform sensitive measurements even 
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in the simplest case of two dotsa. 



Resonant tunneling in arrays of quantum dots and layered semiconductor heterostructures exhibit some similarities. 
The latter situation has been intensively studied in the context of possible Bloch oscillations^. However, Coulomb 
blockade dominates the properties of the arrays of-, Quantum dots so that electron-electron interaction cannot be 
neglected as in the case of layered heterostructuresaja. A way to circumvent this difficulty is to perform an exact 
diagonalization of electron states in the array of coupled dots. Then one considers independent tunneling transitions 
OS ] between the resulting many-electron stateso. This we call the independent channel approximation. This is approximate 
because it disregards the simultaneous tunneling an electron through multiple conduction channels. Another approach 
is to restrict the basis to the resonant states of the uncoupled dots. Then the tunneling between the dots and the 
reservoirs is incorporated ffito a modified Liouville equation for the density matrix in this basis. For two quantum 
dots this has been done inQ and here we extend this approach to the case of an array of an arbitrary number of dots. 
In this paper we concentrate on an array of quantum dots where dots are connected in series and a tunnel coupling 
■ exists only between neighboring dots. This is the most interesting case because there is a unique path for the current 
and changes in any dot strongly affect the transport through the whole array. The first and last dot of the array 
are connected to leads. We assume that the voltage bias is sufficiently high so that the energy change during the 
tunneling of an electron between a reservoir and the array is much larger than the energy uncertainty due to this 
tunneling. We also assume that the resonant electronic energies in the array lie well between the Fermi levels of the 
leads. This enables us to use the density matrix approach. We consider two limiting cases of the electron-electron 
interaction within the array. In the first case we assume that the long range Coulomb repulsion between electrons 
in different dots of the array is so strong so that only one or no extra electrons are present in the array (Coulomb 
blockade). This is to be contrasted with the case of "free" electrons. As we explain below (Section |l]) we do not 
disregard interactions completely in the latter case but rather account only for strong repulsion within each dot. 

Using the density matrix approach in the basis of localized states we have obtained analytical results for the 
stationary current. Our results hold for arbitrary values of the parameters (within the applicability of our model) 
characterizing the array like dot energies and tunnel couplings: no assumption about homogeneity of the array has 
been made. This may facilitate the comparison with experiments and the design of resonant tunneling devices. 
We report the effects of localization and Coulomb repulsion on the resonant current when the energy level of the 
first and the last dot are independently varied. We have also considered another picture of the transport using the 
approximation of independent conduction channels in the array of dots. Using the density matrix approach in the 
basis of delocalized states we have calculated the occupations of the channels and their contributions to the current. 
We discuss in detail the range of validity of this approximation. In the limit of both weak and strong coupling to one 
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or both of the leads we obtain results in agreement with the former more general calculations. However, there can be 
substantial deviations from the predictions of this model when the tunnel rates and the coherent interdot couplings 
are comparable. To illustrate this we study the dependence of the current on the transparencies of the tunnel barriers 
and find unusual features due to the interference of electrons during tunneling. 

The outline of the paper is as follows: in Section [H] we introduce the density matrix description of a multidot system 
coupled to leads which we apply to the Coulomb blockade case in Section [II and to the "free" electron case in Section 
[V|. In Section we compare the results with those obtained from the independent channel approximation and we 
discuss the deviations. We formulate our conclusions in Section [vj. 



II. ARRAY OF QUANTUM DOTS COUPLED TO RESERVOIRS 

Let's first consider an array of N quantum dots (the "device" ) without any contacts (Fig. [j]).We consider a quantum 
dot as some complex system with discrete many-electron states of which only two ground states participate in resonant 
transport. The one state is related to the other by the addition of an extra electron costing an addition energy which 
includes intradot charging i.e. the-Coulomb interactions between electrons in the dot. (For simplicity we disregard the 
electron's spin degrees of freedomla). We label the resonant states of dot i = 1, . . . , N by the number of extra electrons 
rii = 0,1 and introduce fermionic operators a\,a,i which create resp. annihilate an extra electron in dot i. The 

many-electron eigenstates of the array of uncoupled dots will be denoted by \ni . . .tin) = \{n-k}) — OfcLi ^l™ fc | {0})- 
Let Si denote the energy for adding an electron to dot i and let uij_ji > be the interdot charging energy due to the 
repulsive Coulomb interaction between the pair of extra electrons in dots i ^ j. A coherent interdot coupling with 
matrix element ti — t* accounts for the tunneling of electrons between dots i and i + 1. We obtain a Hubbard- type 
Hamiltonian H = Hq + H u for the array of coupled dots: 

JV N-l 

H = ^ + ^ ti (al+i^ + a\d i+ ^j (la) 

i=l i=l 
N 

H u = ^ ^\i-3\^\o-i^jO-i ( lb ) 
i<j=2 

The density operator of the device a evolves according to the Liouville equation <9 t cr = — i[H, a] (h = 1) which 
conserves probability i.e. Tra(t) — 1. By expanding a in the many-electron eigenstates of the uncoupled dots we 
obtain a 2^ x 2 N density-matrix <x = (o~{n k },{n' k }) : 

o= CT {n k },{n' h }\{nk}){{n' k }\ 
{«fc},{«fc} 

We additionally introduce an iV x N Hermitian matrix p with expectation values of single-electron operators = 
1 V): 

pij = (ajdj) = Tra]ajCT 
Using the fermionic commutation relations we find 

Pii — ^ ' °Vii...l...njv,ni...l...niv (^ a ) 

{n k7 ti} 

Pij = ^ ("I) + 3 1 ^n- A:... n \ .n ...0... I \ (2b) 

where {rik^i.j} indicates that we sum over rik = 0, 1 for all dots fc = 1, ...,7V except We will refer to p as 
the average occupation matrix w.r.t. electrons in individual dots whereas cr is the probability-density matrix w.r.t. 
many-electron states of the array. 

Now we include leads L and R connected to resp. the first and last dot of the device by a tunnel barrier (Fig. 
0). The leads are considered here as electron reservoirs at zero temperature with a continuum of states filled up to 
their resp. Fermi levels pl and pr. We assume that the energy levels of the device are located well between the 
chemical potentials of the leads i.e. pl ~3> £i,U ^> pa (large bias) and that the level widths are much smaller than 
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the bias i.e. — 3> ^l,r (discpte states). Under these assumptions an evolution equation for the density-matrix 
cr of the device can be derivedliBjIj by incorporating the details of the lead states into tunnel rates. Due to the 
high bias electrons only tunnel through the barrier from lead L to dot 1 with a rate Tl = 2-kDl 1 2^1 2 and from dot 
N to lead R with rate Tr = 2itDr \ t_n\ . We assume that the density of states Dl r in leads is constant and that 
the tunneling matrix elements t^.R between lead and dot states depend only weakly on the energy. Due to the 
destructive interference of an electron tunneling between a discrete state in the device and the continuum of states 
in a reservoir the rate for tunneling in .one direction [L — ► 1, N — ► R) is constant in time whereas the rate of the 
reversed process (L <— 1, N <— R) is zercO. In general transitions between discrete many-electron states a, b of a device 
with Hamiltonian H induced by tunneling to and from reservoirs can be included into a modification of the Liouville 
equatiorO: 

d t <Jab = -i[H,a] a b 

— OCT"! Ta^a' + ^ Fb->b')<Tab + ^ ^ ab^-a'b' <^ a 'b' (3) 
a' b' a'b' 

The first term of the modification describes the separate decay of states a — > a' and b — > b' due to (in general) different 
tunneling events between the device and the reservoir with resp. rates T a ^ a / and Tb-*b'- The second term describes 
the joint generation of states a <— a' and b <— b' due to a single tunneling event occurring with rate T a b^ a >b r - the 
coherence lost by the simultaneous decay of states a',b' is transferred to states a, b. When we assume that there is 
a unique path for the current, then each state a is generated from a unique state a' by the tunneling of an electron 
to or from a reservoir: T aa t- a 'b' = ^a^a'0~a'b' ■ The modified Liouville equation conserves probability: summing the 
equations for the diagonal elements gives dt Tr a (t) = so 

We consider the following two cases for resonant transport through the array: (i) "Free" electron (F) case where 
interdot Coulomb repulsion is negligible i.e. Uu-ji = 0; up to N electrons can populate the array and all 2 N many- 
electron states lie between the chemical potentials pi > [ir and participate in resonant transport, (ii) Coulomb 
blockade (CB) case where the interdot Coulomb repulsion is so strong that the smallest charging energy is large 
relative to the bias i.e. un-i S> Ml — A*Rj many-electron states with more than one electron are highly improbable 
and can be neglected for resonant transport. We point out that in general the presence of electrons in the array modifies 
the rates for tunneling to or from the leads by Coulomb repulsionllJ. However, in the limiting cases considered here 
there is either no repulsion (F) or no other electron present in the array (CB) so two parameters Yl and Tr suffice 
to incorporate the details of the electronic states in the resp. leads. The current flowing from dot N to reservoir R is 
determined by the average occupation and the tunnel rate: 

Ij ^- = r RPNN (t) (5) 

e 

With the leads included the dynamics of the average occupation matrix p should be calculated from the full density 
matrix cr which evolves according to an equation of the type ([?]). However, in both cases considered here we can 
derive a dynamical equation for the matrix p which is solved more easily. 



III. CURRENT IN COULOMB BLOCKADE CASE 

When the long range Coulomb repulsion is so strong that at most one electron can be present in the array of dots 
we can restrict the set of many-electron basis states to . . . 0), jO . . . 1 . . . 0),i = 1, . . . , N. The average occupations 

i 

of individual dots are simply equal to the non-zero probability densities of these states and 

Pij — 0o...i...o,o...i...o 

* i 

Conservation of probability (||) suggests that we additionally define an average occupation of the many-electron 
vacuum state: 

Poo = co...o,o...o 
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This quantity is positive < poo < 1 and satisfies a conservation law: 

N 



Poo + Pkk = 1 (6) 



fc=i 

In the restricted basis the matrix elements of the Hamiltonian H which describe the coherent part of the evolution 
of the state only involve Ho (eq. ([la])). Modifying the Liouville equation according to (|^) we obtain a dynamical 
equation for the average occupation matrix 

dtPoo = -r L /5 o + ^rPnn (7a) 
dtPu = i (ti-ipii-i + tipu+i — ti-ip%-u — tiPi+u) 

+r L p 00 Sii - T R p NN 5 iN (7b) 
dtPij = i {sj - pij 

+i (tj-lPij-X + tjPij + l — ti-lPi-lj — tiPi+lj) 

-■^RPiN^jN (7c) 

where j > i = 1,...,N and to = = 0. In the rhs of ( ]7a| ) the negative contribution describes the decay of the 
vacuum state due to the tunneling of an electron from lead L to dot 1 with rate Tl whereas the positive contribution 
describes the generation of this state due to the tunneling of the (only) electron in the device from dot N to lead 
R. In d7c|) there is only a negative contribution due to the tunneling of an electron out of dot N to reservoir R with 
rate Tr. There is no negative contribution with F^. because we have incorporated Coulomb blockade: we disregard 
the decay of many-electron states with 1 electron to a state with 2 electrons which occurs when an electron tunnels 
from lead L into dot 1. We can eliminate poo from eqns. (^) using eq. (^|) and obtain iV 2 equations for the average 
occupations. 

Eqns. (f7j) can be used to describe non-stationary transport with a typical relaxation time scale T^ 1 ^. Here we are 
interested in the stationary limit dtp — only. In general the solution of eq. ^ can be obtained by inverting a 
matrix of dimension ./V 2 . However, since the system under consideration is an array with subsequent tunnel coupling 
most of the equations only relate matrix elements of p on neighboring dots. One can obtain the solution by iterativcly 
expressing all matrix elements in terms of poo and finally making use of the normalization constraint (^). The resulting 
stationary current in general reads as 



1 + 1 F + Tr F 



(8) 



Here dimensionless expression Fn depends only on e±, . . .ejv through the differences £y = £j — Ej,i < j and on 
t\ . . . t N— i . Note the curious property of eq. (||): Fn enters the expressions for both In and In+i- This helps when 
deriving expressions for the current through an array with an increasing number of dots. For N = 2 we reproduce 
the result of Stoof and NazarovEj: 
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(9) 



For N = 3 we obtain for the current through a triple dot system 
7 3 CB _ 1 



I < 1 f 3 | £i 2 £l3 + 2f i - tj £i 2 £l 3 ±tj^j4 ] £ 12 + ^23 £lj + £ 2 3 ^ | Fg ^ , -U 



Tl Tfl V tit2 t\ti ti ti j Ati \ t\ 



(10) 



One can show F^ [ex, . . . g jy, t\ . . . t^r-i) = N for £y <C ti = t. This property will be used later on and it is derived 



in an other way in Section V A . 



We consider eq. (ra) in more detail for equal interdot couplings ti — t and several different configurations of the 



levels £j. At resonance £y = we find 
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and the current is maximal: 



rCB 
-'W, max 



T- L + Y- R N + ^ (N - 1] 



1 

T- L + T- R N + ^ {N - 1] 



(ii) 



Clearly the current is reduced when the array increases in size: the number of states participating in transport 
relative to the number of states of the array decreases oc 1/N due to the Coulomb blockade. In an actual system we 
expect a positive deviation from this decrease to occur when the spatial size of the array exceeds some range over 
which the Coulomb repulsion cannot exclude the occupation of a second dot in the array. Away from resonance i.e. 
Eij = O (e) ^> t the conduction of the device decays exponentially with the size of the array N due to the localization 
of the electron in one of the individual dots: i^p/ e oc T ^(e/t)^ 2 ^^ 1 ^. To illustrate this we vary only the last level 
(Fig. ||a) i.e. we consider the solution of (fjj) for £j = Eat^ 



'JiN- 



1 



J_ + J L( AI + (w _ 1) (£ f )») + I| (W _ 1) 



As we increase the number of dots the curve keeps its Lorentzian shape w.r.t. ejq and always depends on both t and 
F# since the electron is localized only just before tunneling out of the array: 



IT N- 



1 



2t 



tGB 

TV, max 



Now we vary only the first level (Fig. 0b) i.e. we consider the solution of (Jfy for = £i£ji : 

/£ B i 

E ki- 

k=l x 1 




r L r 



For large N the normalized current vanishes when the detuning exceeds the tunnel coupling (since an electron is 
localized in the first dot) whereas near resonance the peak takes on a parabolic shape which is independent of the 
tunnel rate r^: 

f 1- (y) 2 £i/t< 1 
\ e 1 /t>l 



IT 

rCB 

AT, max 



For the case where the energies are configured as a "Stark ladder" of total width e we have plotted the current in Fig. 
|3|b. The localization of electrons clearly dominates the current since the tails of the current peak decrease rapidly 
with increasing N as in Fig. 0a. 



IV. CURRENT IN "FREE" ELECTRON CASE 



When interdot Coulomb repulsion is altogether disregarded all 2 N many-electron states of the array of dots 

must be taken into account. Modifying the Liouville equation with H = Hq according to the general prescription (|^) 
we obtain the following set of 2 2N equations for the density matrix: 



dtO- ni . 



= -i H ,<j 



^ r i((l - n l) + (1 - «-i)K ll ...n«,n' 1 . 



+ T R (l-n N )(l 

) . . .n jv - i 1 ,7i^ . . .n' N _ 1 1 



(12a) 

(12b) 
(12c) 
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Here (12a) describes transitions between the non-orthogonal states of the device with a fixed number of electrons 
due to the tunneling between neighboring dots. In co ntras t to the Coulomb blockade case both tunneling into (T^) 
and out of the array (I\r) give negative contributions (|l2bj). Furthermore, there are tunn el processes which induce a 



transition between 2 pairs of many-electron states and give a positive contribution (12c) to the coherences. Only a 
subset of (|l^) forms a closed system of equations for the diagonal and some non-diagonal elements of er (The remaining 
equations only couple a closed set of non-diagonal elements which are irrelevant). From this subset we can derive 
an even simpler dynamical equation for the average occupation matrix p. Let's first derive the coherent part of this 
equation by considering "free" electrons in the array of dots without the leads. Because the Hamiltonian consist of 
only one-electron operators (Hq — Ylkl=i ^okia^ai) , the commutator in the Heisenberg eq. of motion for pij — (ajai) 
is readily expressed in other one-electron operators and the matrix elements of Hq: 

idtptj = (\a]oi,H ) 

N N 

= £ H 0il (flXai) - J2(al ai )H okj = [H , p] y (13) 



1=1 k=l 



wher e i, j — 1,...,N. This equation can also be derived by taking the average (defined in eq. (g)) of the coherent part 



(12a). The contributions which describe the coupling to the leads are found by adding the average of the incoherent 
contributions ( |l2b| ) and (12c) to the rhs of (|l3|). We obtain the following closed system of only iV 2 equations which 



describes the dynamics of the average occupation number matrix: 

dtPii — i(ti-lPii-l + tiPii+l — tiPi+li ~ h-lpi-li) 

+T L (1 - pn)Sn - T R p NN S lN (14a) 
d t pi 3 = i {sj - Si) pij 

+i (tj-ipij-i + tjPij+i — ti-ipi-ij — tipi+ij) 

--j^LPijSn - ^ R p iN S jN (14b) 

where j > i = 1, ...,N and to = = 0. These equations closely resemble those for the Coulomb blockade case: 
t he c oherent contributions are exactly the same, but in contrast to eq. ( ffh] ) the average occ upat ions Ph,pnn in eq. 



(14a) are not coupled by incoherent transitions to some vacuum state. Furthermore, in eq. (14b) there is a negative 
contribution due to the tunneling of an electron into dot 1 with rate Tl which is absent in eq. 
blockade. 



dUIi 

(7q) due to Coulomb 



Despite the close resemblance to eq. (|7|) the (calculation of) stationary solution of eqns. (14) for the general case 
is far more complicated. An analytical expression seems neither feasible nor instructive and we consider only a few 
representative cases here. Assuming equal couplings U — t one finds that the resonant current peak (ey = 0) is 
independent of N: 



T F 1 

1 AT.ma.if J- 



(15) 



e J_ j_ r, 

Tl Tr 4t 2 ^ 4f, 2 

This result was previously obtained by Frishman and Gurvitz for tunneling through multifile-well hctcrostructuresi 
using essentially the same approach as for the derivation of the modified Liouville nmaliorJlil which we have used (eq. 
(0)). By solving eqns. (|TJ) for N — 2 we reproduce the result obtained by GurvitzE3,E3 from eqns. ( fO ) 
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(16) 



where r = Tl + Tr. For N — 3 solving the 9 eqns. (|Ij) gives the current through a triple dot: 

1 = 1 

e / \ 2 / \ 2 / \ 2 I -L.L12. 3_*l£23. _ 1_ ( *i _ *2 | £13 

1 , 1 1 / £11 \ 1 / £23 \ l/^i h \ Tl_ Tr \r L ti r K t 2 t 2 r \t 2 ti J t 2 

Tl r R r L \t! J rR\t 2 ) r\t 2 tj \t\ \t\ j_ ^_ + i jz 



r L ^ r R \ t 2 1 ^ r V *2 I 2tf 

(17) 
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If one would first calculate <r then the subset of 20 relevant equations of ( |l2| ) (containing in total 64 equations) needs 
to be solved. For the case where the energies are configured as a "Stark ladder" of total width e we have plotted the 
current in Fig. ||a. Comparison with Fig. [^) shows that the current for the "free" electron case is less sensitive to 
localization effects than in the Coulomb blockade case. 



V. APPROXIMATION OF INDEPENDENT CONDUCTION CHANNELS 



It is instructive to consider an approximate approach of independent conduction channels to our problem. In Section 



V A this approximation is introduced assuming that both reservoirs are weakly coupled to dots in the array. In this 
case the time needed for tunneling to or from a reservoir T^ 1 is much longer than typical time t -1 of the evolution of a 
coherent state in the array. Therefore an electron completes many coherent oscillations in the array before tunneling. 
Somewhat surprisingly this approach ca n als o be used in the opposite limit of very strong coupling to the reservoirs 
as is shown in Section VB. In Section VC we compare the results of the independent channel approximation and 
discuss some peculiar features of the more general results obtained in III and |lV|. 



A. Weak coupling to the leads 

We first consider the array of dots without the leads. Let's denote the N localized (delocalized) eigenstates of a 
single electron in the array of uncoupled (coupled) dots by \i) where i = 1, . . . , N and the single-electron vacuum 
by |0). By transforming to the basis of delocalized states that diagonalizes Hq (eq. (p"aj)) we obtain new fermionic 
operators 

N N 

a? = EW4> a* = EW 3 * 

k=l k=l 

Since in a delocalized state there is a non-zero probability for finding an electron in both dot 1 and dot N such a state 
can be regarded as a conduction channel which carries a current. The new operators add an electron to channel 
resp. remove an electron from channel \i'). By expanding the density operator of the array a in the "many-channel" 
basis states \n\ . . . n^') = = Y\k=i Q'k™'' IW we obtain a new density matrix <x'. Like in Section || we define 

an average occupation matrix for the channels p'^ = (a'jd' i ) which can be expressed in the density matrix er'. The 
advantage of the new basis is that the dynamical equations for the average occupations p' ti are decoupled from the 
non-diagonal elements when we only take Hq into account: d t p\j — — \{e\ — £j)p'ij where e^j are the energies of the 

delocalized states i,j. In the Coulomb blockade case the effect of H u is easily translated since the basis transformation 
preserves the trace of p: the total occupancy of the channels is restricted to values < 1 whereas in the "free" electron 
case all N channels can be occupied (Fig. ||). 

Now we include the coupling to the reservoirs by a Golden Rule approach. Each channel \i') is connected to lead 
L with modified matrix element t\ = (i'\l)tL and to lead R with t l R — (i'\N)tn: an electron in lead L can tunnel 
through any channel to lead R. The rate for tunneling into channel is proportional to the probability to find 
the electron in dot lwhereas the rate for tunneling out of channel is proportional to the probability to find the 
electron in N: 

ri = r L |(z'|i)| 2 r^ = r«|<z'|A0| 2 (is) 

These probabilities depend on the energies £j of the localized states and the matrix elements tj which couple them. 
When we assume that the uncertainty in the energy during the tunneling of an electron between leads and dots is 
much smaller than the level splitting we can disregard correlations due to the simultaneous tunneling of one electron 
through multiple channels. Thus for weak coupling we can include the leads by only modifying the equations for 
the diagonal elements of p' i.e. the channel occupations. The channels can be treated as independent and their 
contributions to the current can simply be added. 

In the Coulomb blockade case (Fig. |5|a) the interdot repulsion couples the occupations of the channels and restricts 
their sum X^=i p'u < 1 • The tunneling of an electron into the empty device with rate T Z L fills a channel \i') and must 
tunnel out with rate T l R before another electron can tunnels into the device and occupy one channel. The equations 
for the occupations are 
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N N 

00 



9 t Poo=E r ^-(E r i>c 

i=l i=l 

ftp-, = rU - iV« 



Using the implied conservation of probability (J2iLo p'u (*) = 1) tne stationary solution (lim t _ 00 9 t/ o^ (t) = 0) is easily 
found. Summation of the contributions of the independent channels gives the stationary current: 



N 



'-f- = E r ^ (oo) = -^g- = l \ (20) 

This is just the expression for the current through one double barrier: the tunneling time is just the sum of the 
tunneling times of the individual barriers where the time for tunneling through the right barrier is increased by a 
factor 

n 2 El(*'|i}| 2 

F N= ^\VM_ > 2=1 = 1 (21) 

i=i \{i'\ N )\ £ |(z'|A0| 2 

i=l 

In the "free" electron case (Fig. ||b) the occupations of the channels are not coupled and their sum is restricted 
only by £ i=1 p'u < N . The device is equivalent to N independent double barriers in parallel where the occupancy 
of channel i obeys 

9tPu = ri (i - P ' u ) - iv« 

From the stationary solution (lim^oo dtp'^ (t) = 0) the current is readily found: 

i=l i=l — ^ H - 

L 1 7?. 

The tunneling time through each double barrier is the sum of the tunneling times (r^ fl ) 1 > r^ 1 ^ whereas the 
tunnel rate of the device is the sum of the tunnel rates of the N double barriers. 

For a few representative cases we have explicitly worked out the approximate approach in the limit of weak coupling 
to both leads Tl R <c £y, h. For the case of N = 2 dots we have first calculated the 2 exact eigenstates of one electron 
in the array of coupled dots without the reservoirs. From these we obtained the tunnel rates to and from each of 
the delocalized states. Finally we summed the contributions of the independent channels to the total current: eq. 
( pp| ) precisely gives the result ^ without the term (rs/2ii) 2 whereas ( p2|) gives the result ( |l6|) without the term 
TiXr/ (2ii) 2 (The latter result was previously found iru for the case of a double quantum well). For the case of TV 
dots with equal interdot couplings tj = t and aligned levels e, = £j a simple result is also possible because in each 
delocalized state the electronic densities in the dots are spatially symmetric. Then the modification factors of the 
rates are equal for each channel (|(«'|1)| 2 = |(i'|A r )| 2 ) and they cancel in ( pp| ) and (p2|). Thus for weak coupling to 
both reservoirs we obtain 

rCB 1 
-"TV. max 1 



(23) 



I F 1 



which is just ( pi] ) and ( p"5| ) without the terms (N — 1) (r^/2t) 2 and T^Tr/ (2i) 2 respectively. Coulomb blockade 
increases the effective time for tunneling out of the device by a factor N w.r.t. the "free" electron case. The current 
increases linearly with each rate Tl,Tr ; <C t as expected from the enhanced tunneling. 
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B. Strong coupling to the leads 



The terms which were found missing above become important when one of the tunnel rates to the reservoirs is 
comparable to the coherent interdot coupling. The correlations between the conduction channels can no longer be 
disregarded in this case: eqns. (Q) resp. ( |l4| ) do not decouple into separate sets of equations for diagonal and non- 
diagonal elements on the basis of delocalized states. However, we will now show that in the limit of strong coupling 
to one or both of the reservoirs we can use the independent channel approximation again. (Note that the energy 
uncertainties are assumed to be smaller than the bias <C hl — Hr) 

First we consider the case where only the last dot N is strongly coupled to the right lead: Tl <C £ij,ti <C T R . The 
electronic state of the dot N and reservoir R form a continuum of states with a Lorentzian spectral density which is 
approximately constant over the energy range t^-i <C T R : 



D n +r (e) = 



1 4 



2?r ( e - e N ) 2 + (T R /2y 2itT r 



The array of the remaining N — 1 dots is weakly coupled to this continuum of states with matrix element i/v-i- 
Therefore we can apply the independent channel approach: the tunnel rate from dot N — 1 to the continuum on the 
right is found with the Golden Rule 

At 2 

f R = 2irt 2 N _ 1 D N+R (e) = — (25) 



For the case of aligned levels = and equal couplings U=t substitution of T R — > T R and N —> N — 1 in (_3|) resp. 
( pi| ) gives the maximum current in the limit of strong coupling to the right lead: 

rCB 1 

^ = ^4^7 = ^-1^ : ^ 

(27) 



4m 1 1 



e J_.J_J_.Efl 
Tl f R T L + At 2 

which is in agreement with our result jlll) resp. (^). For strong coupling to the right lead the current decreases as 
r^, 1 which is somewhat surprising because tunneling is enhanced. The origin of this effect is the formation of linear 
combinations of the discrete state in dot iV with reservoir states with energies roughly Jpetween en ± Fr- Because the 
tunnel processes from reservoir states back into the discrete state destructively interfereO, the discrete state irreversibly 
decays into the continuum. The resulting spectral density decreases with the energy uncertainty T R . The eventual 
decrease of the-currcnt with the tunnel rate T R can also be interpreted as a manifestation of the quantum Zeno effect 
as discussed irc_. More general we have the following relation: for <c £y , U <SC T R 



2 



I N (ej . . . en, ti . . . ijv-i; r_, T R ) = 7 jv _ 1 (ei . . .£jv-ij*i • • -tN-2'^L, — — 
ijj [e\ . . . en, t\... ijv-ij Ti, T R ) = iJ i _ 1 (ei . . . gj^-ij*! • ■ ■ *JV-2i Tl, 



At 2 



This is clearly satisfied by ([lO]) resp. and (§/_) and the general form (||) has this property. 

Now consider the case where the array is coupled strongly to both leads i.e. ey,tj <C r^^. In the "free" electron 
case the reservoir L coupled to dot 1 gives a new continuum of states with a spectral density which is approximately 
constant over an energy range t\ <C Tl- 

1+L{£) 2 H£ - £l) 2 + {TL /2f~^r L 

The tunnel rate from this continuum to dot 2 is 

At 2 

T L = 2irt 2 D 1+L (e) = -± 

1 L 







The remaining N — 2 dots are weakly coupled to a continuum on the left with matrix element t\ and to the right 
with ijv-i and we can applying the independent channel approach to the N — 2 conduction channels. For the case 
of aligned levels £j = £j and equal couplings ti — t substitution of both Fl — > fx, and — > Tr in (24) gives the 
maximum current in the limit of strong coupling to both leads 



r l l 

max 



i i r, 



t l r R 4i 2 

in agreement with result (^5[). More general we have the relation for £j,tj <C Tl r: 



^2 4^2^ 

Ijv (ei . . . en, t\ . . . tjsr-VjTL, Tr) = I N _ 2 { £ 2 ■ ■ ■ £jv-2, • • • £jv-2; ■=r- ! — p — , 



r 



L 



which is satisfied by (|T^)- In the Coulomb blockade case interdot repulsion prevents the discrete state in dot 1 from 
mixing with the reservoir: even if the energy uncertainty allows tunneling into dot 1 to occur on a small time scale 
r^ 1 , the next electron will have to wait for the previous one to tunnel out which occurs on the much larger time scale 
max{e t j 1 , t^ 1 , T^ 1 }. In the limit of strong coupling to both leads the current in this case is independent of Tl and is 
correctly given by ( p6| ) in agreement with result dTT|). 

Finally, we consider the case where only dot 1 is strongly coupled to the left lead i.e. <§; £ij,U <C Tl- In 
the Coulomb blockade case the weak coupling result still applies as explained above. In the "free" electron case the 
discussion is completely analogous to the case of strong coupling of dot N to the right lead and the result is obtained 
by simply interchanging L <-> R and 1 <-+ N.. 



C. Intermediate coupling to the leads; maximum current 



The competition between enhanced tunneling for weak coupling to the leads and destructive interference in the 
opposite limit implies that the current reaches a maximum value when the rate for tunneling into a dot is comparable 
to the coherent coupling to th e ne ighboring dot. We can find the precise location and value of this maximum with 



the results obtained in Section III and IV which also hold in this intermediate case 



First we consider the resonant current peak (Jill) and (15) as a function of the transparency Tr of the right tunnel 



barrier as plotted in Fig. |4|. Starting from zero the current initially increases linearly as expected from the enhanced 
tunneling to the right lead. Then a maximum is reached: 



job l 



r i 

max _ 



r L t 



r * = VA^l 2t ^ 
T R = 2t (28b) 



Increasing the transparency further will reduce the current as explained above. The maximum occurs when there is 
an optimal balance of the coherent tunneling into dot N and incoherent tunneling from this dot to reservoir R. At 
this point the effective time for tunneling out of the device is the same for weak (eq. ( |23| ) resp. (p4|)) and strong 
coupling to the right lead (eq. ( |26| ) resp. (p7|)) i.e. at these values of Tr we have 

- 1) (CB) 

(F) 

where Tr is given by (|25|). For the case of "free" electrons the effective time for tunneling out and therefore the 
position of the maximum is independent of the number of dctts. N. (For this case the Tr dependence of the resonant 
current peak has been discussed for a double dot system inliil) . In the Coulomb blockade case the non-monotonic 
variation of the current with Tr persists for all N. The maximum occurs at a slightly higher value of Tr than for 
the "free" electron case: at Tr = 2t = Tr the effective tunneling time for weak coupling is still larger than for strong 





1 


f R N = 






Fb 


1 


1 




Tr~ 
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coupling because the fraction of channels exclu ded by Cou lomb repulsion is larger for N dots then for N — 1 dots. 
The tunnel rate must be increased by a factor yj N/(N — 1) to exactly balance the effective tunneling times. For large 
N this difference becomes negligible and the maximum occurs at the same position as for the "free" electron case but 
with a much smaller amplitude (Fig. ^). 

Next we consider the resonant current peak ( p"T| ) and (O) as a function of the transparency of the left tunnel 
barrier. In the "free" electron case the maximum current (|lq) remains unchanged when we interchange and Tr. 
Therefore the resonant current peak ([l5]) also displays a maximum as a function of Tl at Tl = 2t. As a function of 
both rates the maximum current is 



iV,max 



t 

2' 



r, 



2t 



As discussed in Section 
resonant current peak (|28e 



Eh there is no such effect in the Coulomb blockade case: due to interdot repulsion the 
will increase with Tl and saturate at a maximal value when Tl^> t 



rCB 

TV, max 



t 



y/N(N-lY 



T L > r fi = 



N 



N - 1 



-2f 



For N = 2 the maximal current as a function of the rates is larger in the Coulomb blockade case whereas for N > 2 
it is larger in the "free" electron case. 



VI. CONCLUSIONS 



We have extended the density matrix approach to resonant tunneling to the case of a linear array of quantum dots 
with strong, long range electron-electron interaction. We have found exact analytical expressions for the stationary 
current in the array for an arbitrary set of parameters (within the applicability of our model) characterizing the array. 
Coulomb repulsion was found to reduce the resonant current by a factor of the order of the number of dots. The 
formation of a localized state in one of the dots when the energy level is displaced results in an exponential decay of the 
current with increasing size of the array. Our approach takes into account correlations between conduction channels 
in the array due to the coupling to the electron reservoirs. This makes our results also valid for relatively strong 
tunnel coupling to the reservoirs where the independent channel approximation does not work. These correlations 
manifests themselves in the eventual decrease of the resonant current when the rate for tunneling into the reservoir is 
increased. 
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FIG. 1. Linear array of TV quantum dots coupled to leads L and R. The energy levels of the uncoupled dots are given by 
full lines when relevant for resonant transport and dashed lines when irrelevant. 

FIG. 2. Normalized resonant current for the Coulomb blockade case with TV = 2, 3, 4, 5, 6 (full curves downwards) and TV — oo 
(thick full curve), (a) Variation of the last last level i.e. e< = £jv<5ijv . Increasing TV does not alter the Lorentzian form of the 
curve, (b) Variation of the first level i.e. Si = EiSu. Electrons are localized in the first dot after tunneling through the left 
barrier resulting in the exponential decay of the current tails with TV. For TV = 2 the curves in (a) and (b) coincide. 

FIG. 3. Normalized current through an array of TV = 2, 3, 4, 5, 6 dots with energies configured as a Stark ladder 
Ei = i/ (TV — l)e of varying width e. (a) "Free" electron case: the curves for TV = 2, 3 (indicated by the arrow) coincide 
for an isotropic array U = t with r_L = Tr. (b) Coulomb blockade case. 

FIG. 4. Maximum resonant current as a function of the transparency of the right barrier: Coulomb blockade case for 
TV = 1, 2, 3, 4, 5, 6 (solid lines downwards) and "free" electron case (dot-dashed line for any TV > 1). The dotted lines show the 
position of the maxima. 

FIG. 5. Tunneling through TV independent channels, (a) Coulomb blockade case: an electron tunneling through one of the 
channels blocks the remaining TV — 1 channels. The state must first decay to the vacuum before another electron can enter one 
of the channels, (b) "Free" electron case: TV parallel channels for tunneling are available. 
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